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ABSTRACT 



Aims. Detect tertiary components of close binaries from spectroscopy and light curve modelling; investigate light-travel time effect 
and the possibility of magnetic activity cycles; measure mass-ratios for unstudied systems and derive absolute parameters. 
Methods. We carried out new photometric and spectroscopic observations of five bright ({V) < 10.5 mag) close eclipsing binaries, 
predominantly in the southern skies. We obtained full Johnson BV light curves, which were modelled with the Wilson-Devinney code. 
Radial velocities were measured with the cross-correlation method using IAU radial velocity standards as spectral templates. Period 
changes were studied with the O-C method, utilising published epochs of minimum light (XY Leo) and ASAS photometry (VZ Lib). 
Results. For three objects (DX Tuc, QY Hya, V870 Ara), absolute parameters have been determined for the first time. We detect 
spectroscopically the tertiary components in XY Leo, VZ Lib and discover one in Q Y Hya. For XY Leo we update the light-time effect 
parameters and detect a secondary periodicity of about 5100 d in the O-C diagram that may hint about the existence of short-period 
magnetic cycles. A combination of recent photometric data shows that the orbital period of the tertiary star in VZ Lib is likely to be 
over 1500 d. QY Hya is a semi-detached X-ray active binary in a triple system with K and M-type components, while V870 Ara is a 
contact binary with the third smallest spectroscopic mass-ratio for a W UMa star to date (q = 0.082 ± 0.030). This small mass-ratio, 
being close to the theoretical minimum for contact binaries, suggests that V870 Ara has the potential of constraining evolutionary 
scenarios of binary mergers. The inferred distances to these systems are compatible with the Hipparcos parallaxes. 

Key words, stars: binaries: close - stars: binaries: eclipsing - stars: binaries: spectroscopic - stars: binaries: general 



1. Introduction 



Modelling variations of close eclipsing binaries is a power- 
ful method for determining fundamental physical parameters of 
stars. This is because the observed brightness, colour and ra- 
dial velocity changes give strong constraints on the geometric 
configuration of a given system. W UMa-type stars are low- 
mass eclipsing contact binaries with orbital periods less than 1 
day, showing continuous light variations. The components fill 
their Roche lobes so that the strongly distorted stars touch each 
other at the inner Lagrangian point. These systems have long 
been known as the most numerous of all stars, with roughly one 
W UMa binary per 500 ordinary dwarf stars (Rucinski 2002a), 
which is one of the reasons why many of them, including bright 
southern ones, are left unstudied until now. In particular, mea- 
suring radial velocities has been a tedious task with very slow 
progress. For instance, Bilir et al. (2005) noted that among the 
751 recorded W UMa binaries in the recently revised fourth 
edition of the General Catalog of Variable Stars, only 129 sys- 
tems were found to have radial velocities, mostly northern ones. 
Very recently, Rucinski & Duerbeck (2006) and Duerbeck & 
Rucinski (2007) published radial velocities for 23 predominantly 
southern close binaries, but apart from that, there is basically no 
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dedicated program for obtaining new data in the south. There 
are several important problems related to the formation, inter- 
nal structure and evolution of these systems, like the kinematics 
(e.g. Bilir et al. 2005), energy transfer between the components 
(Csizmadia & Klagyivik 2004 and references therein) magnetic 
activity and its cyclic nature (Borkovits et al. 2005) and the fre- 
quency of additional components of W UMa stars (Pribulla & 
Rucinski 2006 and references therein), all in need for more ex- 
tensive data. The latter one is particularly important, because it 
is the angular momentum transfer in hierarchic triple systems 
that can lead to such close binary systems as W UMa-type stars, 
so one can even hypothesize that all contact binaries reside in 
triple or multiple systems (Pribulla & Rucinski 2006; D' Angelo, 
van Kerkwijk & Rucinski 2006). Testing this hypothesis is im- 
portant, because angular momentum evolution, in which mag- 
netic braking was shown to be a key factor to explain the ob- 
served properties of contact binaries (Mochnacki 1981; Stepien 
1995), can be strongly affected by the presence of a third body 
(Eggleton & Kiseleva-Eggleton 2001). In a wide range of initial 
conditions an originally detached binary can reach the contact 
phase and this process can be altered by a tertiary star in the 
system (for a recent review see Eggleton 2006). 

Late-type contact binary stars are known to be very ac- 
tive objects with stellar spots (Maceroni & van't Veer 1996) 
and X-ray radiation due to chromospheric activity (Stepien et 
al. 2001; Chen et al. 2006; Geske et al. 2006). These pieces 
of evidence can be interpreted by the Applegate-mechanism 
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Table 1. New times of primary (I) and secondary (II) minima 
for the target stars. 



Trrm [HJD] 


r mm [HJD] 


r raln [HJD] 




XY Leo 




2453412.1839 II 


2453415.0236 II 


2453416.1621 II 


2453413.0367 II 


2453415.1671 I 


2453435.0538 I 


2453413.1793 I 


2453416.0209 I 


2453436.0479 II 




VZLib 




2453189.0102 II 


2453190.9776 I 






DXTuc 




2453187.1669 II 


2453188.2989 II 


2453189.2421 I 




QYHya 




2453186.9081 I 


2453187.9276 II 






V870 Ara 




2453185.1379 I 


2453195.1314 I 


2453196.1325 II 
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(Applegate 1992), which involves the interchange of magnetic 
and kinetic energy. An observable manifestation could be cyclic 
orbital period variation, which seems to be present in many con- 
tact binaries, although not necessarily caused by magnetic activ- 
ity - for instance, light-travel time effect due to a third body can 
also cause cyclic period change. Hence one must take caution 
when interpreting the O-C diagrams in terms of cyclic changes. 
Observations so far seem to agree with the predictions of the 
orbital period-magnetic modulation period relation (Lanza & 
Rodono 1999) and some of the studied systems were also anal- 
ysed from this point of view. 

Here we report on multicolour and radial velocity measure- 
ments of five close eclipsing binaries, of which three have never 
been studied before. In addition, the sample includes the well- 
studied quadruple system XY Leonis, which is one of the most 
fascinating cases of light-travel time effect in a periodic variable 
star (Gehlich, Prolls & Wehmeyer 1972; Yakut et al. 2003), and 
the triple system VZ Lib. Besides revealing new information on 
these two objects, the observations also allowed a comparison of 
our results with previous studies. 

2. Observations and analysis 

The observations were carried out in three observing runs at the 
Siding Spring Observatory, Australia. We took two-colour pho- 
tometry and optical spectroscopy on six consecutive nights be- 
tween June 28 and July 4, 2004 for V870 Ara, QY Hya, VZ Lib 
and DX Tuc. Data for XY Leo were gathered on 4 nights in 
February, 2004 and 7 nights in March, 2005. The photometry 
was done using the 40-inch telescope of the Australian National 
University (ANU) in Siding Spring, equipped with the Imager 
CCD detector (2148x2048 pixels) and Johnson B and V filters. 
The exposure times ranged from 3 to 15 seconds, depending 
on the target brightness and weather conditions. For the optical 
spectroscopy, we used the ANU 2.3-m telescope and the Double 
Beam Spectrograph, recording the second order spectra of the 
1200 mtrr 1 grating, which gave a nominal spectral resolution of 
/I/ A/1 as 7000 at the Ha line. We used only the red beam, cover- 
ing ~ 1000 A between the sodium D doublet and the Ha line. 

All data were reduced in a standard fashion using IRAH3- 
Photometric reductions included bias and fiat-field correction, 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Fig. 1. A simulated CCF profile and the three-component 
Gaussian fit (see text for details). 



while instrumental magnitudes were determined with aperture 
photometry, because the stellar fields around the programme 
stars were quite empty, so that accurate psf-fitting was not pos- 
sible. The instrumental data were tied to the standard system by 
observing equatorial photometric standard stars from the list of 
Landolt (1992), except for XY Leo, which was only observed 
under non-photometric conditions. The estimated uncertainty of 
the absolute values of the B and V magnitudes is about +0.01 
mag. 

From the light curves we determined epochs of minimum 
light by fitting low-order polynomials. We measured 19 new 
times of minimum, which are collected in Table Q] The typical 
uncertainty is about +3 x 10~ 4 d. With these we updated the O- 
C diagrams that were based on all published observations in the 
literature. For V870 Ara, QY Hya and DX Tuc we could derive 
a more accurate orbital period. 

The spectra were also reduced with standard IRAF tasks, 
including bias and fiat-field corrections, cosmic ray removal, 
extraction of one-dimensional spectra, wavelength calibration 
and continuum normalization. The exposure times ranged from 
2 to 5 minutes, and NeAr spectral lamp exposures were regu- 
larly taken to monitor the wavelength shifts of the CCD spec- 
tra. Radial velocities were determined by cross-correlation, with 
IAU radial velocity standards HD 187691 (F8V, V r = -0.2 
km s 1 ) and HD 80170 (K2III, V r = 0.5 km s 1 ). The calcu- 
lated cross-correlation functions (CCFs) were fitted with two- or 
three-component Gaussians (in case of detected tertiary compo- 
nents), whose centroids gave the radial velocities for each com- 
ponents. In most cases the Gaussian approximation gave reason- 
ably good fits for the central parts of the CCFs, while departures 
from the Gaussian shape typically occured only 200-300 km s _1 
away from the maxima. 

The shape of the CCF clearly showed the presence of a ter- 
tiary component in three of the target stars. In those cases we 
recorded its mean radial velocity, which can be used in the fu- 
ture to detect orbital motion around the barycentre of the triple 
system. The binary radial velocity curves were fitted with sine 
curves, which gave the mean velocity of the binary (V y ), the or- 
bital velocities and the mass-ratio of the components (K\, K2, 
q). These were then used as fixed input parameters for the light 
curve models. 

The rms scatter of the derived velocity data around sine- 
wave fits suggests that the typical velocity precision is about +10 
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km s _1 per point. Although Rucinski (2002b) demonstrated that 
the use of the broadening function (BF) for determination of ra- 
dial velocities of contact binaries is superior to the CCF method, 
we chose to apply the CCF for two reasons. Firstly, IDL was 
not available to us and thus we could not use the BF codefl 
Secondly, the data have relatively poor spectral resolution and 
thus the CCF approach is sufficient. Nevertheless, besides esti- 
mating probable errors from the data, we have also tested the 
accuracy of the velocity measurement method in the following 
way. First, a theoretical spectrum of a T e ff = 5900 K, log g = 4.0 
star was calculated with Robert Kurucz's ATLAS9 code between 
5800 A and 6800 A, and resampled to a resolution of ~ 8000. 
Then, a contact binary broadening profile at quadrature phase 
(</> = 0.25) was generated by Slavek Rucinski's WUMA4 code, 
applying the physical parameters of VZ Lib, one of our program 
stars (see Sect. 3.2). The model spectrum was then convolved 
with the broadening profile to get a contact binary model spec- 
trum. A scaled original model spectrum (scaled to 20% of the 
initial flux level) at zero velocity was added to the contact binary 
model spectrum to mimic the presence of the tertiary component. 
This binary+tertiary spectrum was then cross-correlated with the 
original model spectrum. The resulting CCF, which is very sim- 
ilar to most of the real CCFs of our program stars, is shown in 

Fig.m 

The radial velocities were determined via fitting three- 
component Gaussians, as described above. The centroids recov- 
ered for the model spectrum are a\ = -80.5 km s _1 , «2 = 259.1 
km s -1 and «3 = 0.0 km s _1 , while the model velocities were 
initially vi = -84 km s _1 , V2 = 252 km s _1 and V3 = 0.0 km s _1 . 
It is seen that the error of the velocity determination is within 
10 km s for the binary components and less than that (proba- 
bly a few km s~') for the sharp-lined third component. These are 
clearly larger than the 1-2 km s _1 error determined by Rucinski 
(2002b) for the BF-method, but our spectra are of lower resolu- 
tion than those applied by Rucinski (2002b). The numbers also 
show that the recovered velocities are shifted to the input val- 
ues in the same direction, which means the effect is systematic, 
affecting both velocity amplitudes. We note that Rucinski and 
co-workers found that CCF tends the reduce the measured ra- 
dial velocity amplitudes relative to the BF results (e.g. Rucinski 
2004). Conservatively, we adopt a ±8 km s _1 error in the derived 
amplitudes K\ and K2 and +2 km s _1 per point as the measure- 
ment uncertainty of the radial velocity of the tertiary component. 
These errors were added in quadrature to the formal errors of the 
calculated fits. 

The light curves were modelled using the 2003 versiorQ of 
the Wilson-Devinney (WD) code (Wilson & Devinney 1971, 
Wilson 1994, Wilson & van Hamme 2003). Besides the spectro- 
scopic parameters, we determined the effective temperature, T\, 
of the component eclipsed in zero phase from the (B - V) - T e ff 
calibration of Gray (1992). The lower limit of the semimajor 
axis, A, was calculated from K\, K2 and P. These input parame- 
ters were kept fixed while running the WD code on the original 
individual data points. Limb-darkening coefficients were taken 
from Diaz-Cordoves, Claret & Gimenez (1995), while the bolo- 
metric albedos and gravity-darkening coefficients were set to 0.5 
(Rucinski 1969) and 0.32 (Lucy 1967), respectively. The third 
light was included in the fitting procedure, except for DX Tuc 
and V870 Ara, for which we do not find evidence for a third 
component. In three cases, the brightness difference of the con- 
secutive maxima (the O'Connell-effect, hereafter denoted by AV 



in the V band) suggested the presence of spots, which were 
added to the solution. The astrographic coordinates of the center 
of a given spot were fixed after several trials. 

In case of the multiple systems, we also estimated the phys- 
ical parameters of the third bodies using the method described 
in Borkovits et al. (2002) and Csizmadia (2005). The light curve 
solutions yield the third light in different colours and the colour 
index of the tertiary component follows from the colour index of 
the primary star as 



(B-V) 3 = (B-V)i-2.5 



logCy^-logC-^) 

L 3,V M,y 



(1) 



2 http://www.astro.utoronto.ca/~rucinski/SVDcookbook.html 

3 ftp://ftp.astro.ufl.edu/pub/wilson 



where L3 « 4nl^ (Z3 is the third light) and the luminosities are 
given in arbitrary units calculated by the WD-code. The colour 
index of the primary was estimated from the determined effective 
temperature using the temperature-(fi - V) relations by Flower 
(1996). With only two colours we could not identify the location 
of the tertiary star in a two-colour diagram. Therefore, we as- 
sumed that the tertiary component is a main-sequence star. Then 
its spectral type can be estimated from the colour index using the 
Bessell (1990) tables, yielding an approximate mass estimate. 

All the standardized light curves and radial velocity 
data are available electronically at the CDS via anony- 
mous ftp to cdsarc.u-strasbg . fr (130. 79. 128. 5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/AA/. . ./. 



3. Results 

3.1. XYLeonis 

XY Leo is a well-studied contact binary which has long 
been known for its cyclic period variations (Gehlich, Prolls & 
Wehmeyer 1972; Yakut et al. 2003). In stark contrast to the mul- 
titude of photometric studies (Yakut et al. 2003 and references 
therein), only two spectroscopic studies were aimed at measur- 
ing radial velocity curves (Hrivnak et al. 1984; Barden 1987). 
XY Leo has four components. One pair is a W-type contact bi- 
nary (i.e. the smaller star is the hotter one) with very short orbital 
period (P * 0.28 d) and late spectral type (Kn). It is a chro- 
mospherically active variable (Vilhu & Rucinski 1985), but the 
activity is largely dominated by the BY Dra-type binary compo- 
nent that consists of two red dwarf stars (Barden 1987). 

Spectra taken in early 2004 showed a similar, though slightly 
weaker and unresolved Ha emission than the one depicted in fig. 
1 of Barden (1987). The cross-correlation profiles of the spectra 
clearly indicated the presence of a third component (left panel 
in Fig. [2j>, so we fitted three Gaussians to measure radial veloci- 
ties of the eclipsing pair (right panel in Fig.|2|. The determined 
spectroscopic elements (Tabled agree within the errorbars with 
the corresponding parameters of Barden (1987), except the y- 
velocity. Since our observations were taken 19 years later than 
those of Barden (1987), which is close to the orbital period of the 
wide binary, this shift in V y might be due to the strong third light 
that distorts the whole CCF profile. Nevertheless, the agreement 
in the measured amplitudes suggests that the resulting velocity 
shift caused by the distorted CCF is very similar for both com- 
ponents. 

The light curves (Fig. [3} do not show O'Connell-effect, thus 
we did not include spots in the light curve fit. In this case we 
had only differential light curves, so the fixed temperature of 
the primary component was estimated from the spectral type. 
The results (Table [2]) are quite similar to those of Yakut et al. 
(2003). The derived colour index of the third component is very 
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Table 2. Physical parameters of XY Leo. Parameters from the literature were taken from Barden (1987 - spectroscopy) and Yakut 
et al. (2003 - light curve fit and absolute parameters). Values in italic denote fixed input parameters. HJD3 and V3 refer to mean 
epoch and radial velocity of the third component. 



Paramptpr 


TfiiQ nanpr 


T ltpratnrp 


Paramptpr 


"This nanpr 


T itpraturp 


Spectroscopy 
A sin 1 [K© J 


1 on 1 c\ 1 1 
l.oU+U. 1 1 


1 Oj4 1 n AQ 

1.84+0.03 


K[ Lkiti/sJ 


11)3. f±yA 


2U4. 1+1.3 


Vy LK1T1/SJ 


— JU. / ±J. L 






1 99 n_i_ 1 n n 




HJD 3 


2453456.5 




V 3 [km/s] 


-70.0+7.5 




Light curve 
1 L J 


fit 

6/. 1+0. 1 


65+1 


111 


A H'trt/^ 1 A AAT7 

4. /206+0.003 / 


4. / i±0.01 


/ 


5.2% 


6.7% 


Q.2 


4.7206 


4.71 


phase shift 


0.0027±0.0002 


— 




0.553+0.009 


0.507+0.043 


1 


1.66(±0.23) 


1.64 




0.522+0.007 


0.483+0.043 


Ti [K] 


4800 


4850 


(/ 3 )j 


0.051+0.005 


0.019+0.007 


r 2 [K] 


4361±8 


4524+14 


a 3 )v 


0.068+0.004 


0.059+0.007 


(B - V)! 


1.02 




(B - V) 3 


2.02 




Absolute 
Mi [Mo] 


parameters 
0.46±0.06 


0.50±0.02 


M r [Mo] 


0.76+0.15 


0.82+0.02 


Ri [Ro] 


0.66±0.10 


0.68+0.02 


#2 [Ro] 


0.83+0.13 


0.85+0.02 


U [L ] 


0.21±0.07 


0.226 


L 2 [L ] 


0.22+0.08 


0.267 


(M bol )! 


6.43±0.23 


6.4 


(M bol ) 2 


6.34+0.23 


6.2 


(M v ), 


6.89±0.23 


6.7 


(M v ) 2 


6.96+0.23 


6.7 


Sp. type (3) 
rf[pc] 


M 

51+6 


631 










v (km/s) Phase 

Fig. 2. Left panel: the CCF-profile of XY Leo in phase 0.75 with the three fitted Gaussians. Right panel: radial velocity curves of 
the eclipsing components of XY Leo. 



red at (B - V)?, = 2.02 mag, being consistent with a BY Dra- 
like M-dwarf binary (Barden 1987). For calculating the distance 
we used the maximum V-brightness of the stars listed in ESA 
(1997). 

With the most extended coverage of the period change, 
we also updated the O-C diagram. Since the publication of 
the Yakut et al. (2003) analysis, new moments of minima 
have been published by Agerer & Hiibscher (2003), Drozdz & 
Ogloza (2005), Hiibscher (2005), Hiibscher, Paschke & Walter 
(2005) and Nelson (2006). We also added a few epochs from 
Tsessevitch (1954a), which were not included in the analysis of 
Yakut et al. (2003) but were important in stretching the time 
coverage as far back as possible. The full O-C diagram of 
XY Leo, plotted in the top panel of Fig. |4] was calculated us- 
ing the ephemeris: HJD min = 2435484.0283 + 0.28410260 x E. 
Then we simultaneously fitted a parabolic and a light-time term 
to calculate the rate of the linear period change and the orbital 
elements of the perturbing binary. For this, we used the same 



code as Ribas, Arenou & Guinan (2002), kindly provided by Dr. 
I. Ribas. 

The results are summarized in Table [3] In general, the 
parabolic term represents a continuous period variation that 
might be explained by some mass-transfer between the com- 
ponents. The period changing rate of XY Leo - 2.85 x 
10~ 8 yr~', which can be translated to a mass-transfer rate of 
-9.6 x 10~ 9 M Q /yr. This period variation rate is in the usual range 
of W UMa systems (see Qian 2001ab; Borkovits et al. 2005 - 
XY Leo was not included in these studies). The secular orbital 
period increase seems to support what Qian (2001a) suggested, 
namely that W-type systems with q < 0.4 and q > 0.4 tend to 
show secular decrease and increase, respectively. 

The light-time parameters in Table [3] give very similar or- 
bital elements for the quadruple system to all the previous inves- 
tigations (Gehlich, Prolls & Wehmeyer 1972; Pan & Cao 1998; 
Yakut et al. 2003), so that we do not re-iterate the discussion of 
the hierarchic component. However, the residuals of the O-C fit 
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Fig. 3. Observed and fitted light curves of XY Leo. 

Table 3. Parameters of the O-C fit. The standard errors in the 
last digits are shown in parentheses. 



Parameter 


value 


Ao-c [d] 


0.02430(1) 


e 


0.11(1) 


6>n 


18.5(5) 


P 3 [yr] 


19.651(2) 


* periastron [HJD] 


2446119(11) 


Pi [d] 


0.28410295(1) 


dP/dt [s/d] 


1.92(6)xl0" 6 



deserve further discussion, because the 201 s rms is larger than 
the typical observational uncertainties. An analysis of the long- 
term period change of the triple W UMa VW Cep revealed sec- 
ondary cycles in the O-C diagram that were interpreted as possi- 
ble evidence for magnetic activity (Kaszas et al. 1998). More re- 
cently, Borkovits et al. (2005) investigated complex period vari- 
ations of five W UMa-type binaries, finding some common fea- 
tures. These include secular period variation at a constant rate 
and a low-amplitude modulation with periods around 18-20 yr 
in four of the five cases. These stars all have spectral type later 
than F8 (so they have convective envelopes), while the fifth one 
is an A-type contact binary. XY Leo is a K-type star, therefore 
the secondary periodicity of the O-C diagram is in agreement 
with the picture suggested by Borkovits et al. (2005) that all 
contact binaries later than F8 have low-amplitude cyclic period 
changes. These can be interpreted as indirect evidence of mag- 
netic cycles in late-type overcontact binaries that are analogous 
to the observed activity cycles in RS CVn systems (Hall 1989). 
The theoretical framework involves the interchange of magnetic 
and kinetic energy (Applegate 1992; Lanza & Rodono 2002), 
predicting time-scales from several years to decades. 

The bottom panel of Fig. [4] shows the best-fit sine wave to 
the residuals of the O-C. Introducing this extra component de- 
creases the rms to ~175 s, which represents a slight but de- 
tectable improvement and is in agreement with the typical pho- 
tometric noise from starspots (Kalimeris, Rovithis-Livaniou & 
Rovithis 2002). The Fourier spectrum of the residuals shows 
that the highest peak has a S/N ratio (Breger et al. 1993) of 3, 
thus its significance is marginal. Nevertheless, the correspond- 
ing modulation period (~5 100 d or 14 yr) and amplitude (0.0016 
d) are very similar to those found by Borkovits et al. (2005) for 
AB And, OO Aql, V566 Oph and U Peg. We therefore conclude 



0.04 




-0.04 1 ' 1 1 1 1 

30000 35000 40000 45000 50000 55000 

HJD (-2400000) 

0.02 I 1 1 1 1 1 



0.01 
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-0.02 1 ' 1 1 1 1 

30000 35000 40000 45000 50000 55000 

HJD (-2400000) 

Fig. 4. Top panel: the O-C diagram of XY Leo with the theoret- 
ical fit. Bottom panel: the residuals and the ~5100 d secondary 
periodicity. 



that XY Leo is likely to bear signatures of short-period magnetic 
cycles in W UMa-type variables, but several years of further con- 
tinuous eclipse timings are needed to improve the significance of 
the detection. 

3.2. VZ Librae 

The variability of the F5-type star VZ Lib was discovered by 
Hoffmeister (1933), who gave neither period nor classification 
but recognized its low-amplitude eclipsing nature. Tsessevich 
(1954b) classified the star as a W UMa-type variable with a pe- 
riod slightly over 8 hours. Claria & Lapasset (1981) reported 
the first photoelectric light curve which showed a difference in 
the eclipse depths for the primary and the secondary minima of 
about 0.1 mag. Interestingly, more recent observations, includ- 
ing ours and those of Zola et al. (2004) find more similar min- 
ima with differences only 0.02-0.03 mag (Zola et al. 2004) and 
<0.01 mag (this paper) in V. There is a relatively bright tertiary 
component in the system discovered by Lu, Rucinski & Ogloza 
(2001). Besides finding radial velocity variations of the third star 
up to 40 km s _1 over a period of 1200 days, they derived a lu- 
minosity ratio of L3/Z42 = 0.20 + 0.04. Contrary to this, Zola 
et al. (2004) arrived to only a few (4-5) percent of third light 
contribution, which was left unexplained. 

The cross-correlation profiles of our spectra clearly indicated 
the presence of a third, narrow-lined companion (left panel in 
Fig. |5]). Its mean radial velocity is -4.8 + 3 km s _1 , which is 
larger than any of the velocities in fig. 5 of Lu, Rucinki & Ogloza 
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Table 4. Physical parameters of VZ Lib. Parameters from the literature were taken from Lu, Rucinski & Ogloza (2001 - spec- 
troscopy) and Zola et al. (2004 - light curve fit and absolute parameters). Values in italic denote fixed input parameters. HJD3 and 
V3 refer to mean epoch and radial velocity of the third component. 
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Literature 
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this paper 
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Ki [km/s] 


84.1+8.2 


68.5+3.8 


V rkmAil 

V y [Kill/ S J 




1+7 ^ 

J 1 . 1 + Z.. J 


K-, Tkm/sl 






HJD 3 


2453189.4 
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fit 
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Z 3 (B) 


0.192+0.007 


0.011+0.005 




0.2777±0.0007 


0.2545+0.0010 


h (V) 
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(B-n 


0.63 




(B - Vh 


0.73 




Absolute 


parameters 
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Mi [Mo] 
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^2 [Ro] 
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U [L ] 
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(Afw)i 
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Fig. 5. Left panel: the CCF-profile of VZ Lib in phase 0.75 with the three fitted Gaussians. Right panel: radial velocity curves of the 
eclipsing components of VZ Lib. 



(2001), ranging from -9 to -50 km s . Apart from pushing 
up the full velocity range of the third component well over 40 
km s" 1 , this value strongly suggests that the 1200 d of obser- 
vations by Lu, Rucinski & Ogloza (2001) was shorter than the 
orbital period of the third component. 

The radial velocity curve (right panel in Fig. |5]l yields a 
higher mass-ratio than that of Lu, Rucinski & Ogloza (2001), 
but it is still compatible with the earlier result of q — 0.24 + 0.07 
that was based on a relatively poor radial velocity solution. The 
light curves (Fig. |6]l are very similar to those of Zola et al. 
(2004), except that the O'Connell-effect was not detectable in 
2004. Consequently, our light curve fit did not include spots on 
any of the components. 

A comparison of the determined parameters with the pub- 
lished ones shows a relatively good agreement (Table 0). The 



differences can largely be traced back to our larger mass-ratio 
and the different amount of the third light. The strong tertiary 
peak in the CCF-profile supports the ~20% luminosity contri- 
bution of the third component, which is also close to the value 
determined by Lu, Rucinski & Ogloza (2001). Since the Zola et 
al. light curve model is not consistent with this, we are confident 
that our light curve solution is more compatible with the avail- 
able spectroscopic information. To remove most of the ambigu- 
ities, a better defined radial velocity curve is highly desirable, 
preferentially from higher-resolution and S/N spectra. 

Because the tertiary component is comparable in luminosity 
(and perhaps in mass) to the eclipsing pair, one can expect mea- 
surable light-time effect in the period change. Since the study of 
Claria & Lapasset (1981) only the Hipparcos photometry, Zola 
et al. (2004) and Krajci (2006) presented new epochs of mini- 
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Table 5. Physical parameters of DX Tuc. 




Fig. 6. Observed and fitted light curves of VZ Lib. 
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Fig. 7. The O-C diagram of VZ Lib between 2000-2006 (see 
text for details). 



mum light, which is hardly enough to reconstruct the hypothetic 
cyclic period changes. On the other hand, the ASAS-3 project 
(Pojmanski 2002) has been observing this system since 2001 and 
produced 259 V-band points scattered between JD 2451918- 
2453277. We combined the recent data as follows. 

First we adopted the linear ephemeris of Zola et al. (2004): 
E = 2452727.4047, P = 0.3582580 d. Then we cut the ASAS 
data into 4 equal subsets, which correspond to four observing 
seasons since the beginning of the project. We phased all the 
subsets using the given ephemeris and the resulting phase dia- 
grams yielded the phase-shifts of the primary minimum. In Fig. 
[7] we plot the results (multiplied by the period), where the last 
two points were calculated from our primary minimum and that 
ofKrajci(2006). 

The shape of the O-C diagram indicates measurable phase 
changes over the last 6 years. This may be interpreted as possi- 
ble evidence for period variations due to light-time effect in the 
triple system, although other mechanisms, like magnetic activity 
or sudden mass exchange via large flares, cannot be ruled out ei- 
ther. Assuming light-time effect, both the long-term radial veloc- 
ity curve (Lu, Rucinski & Ogloza 2001) and Fig. [7] suggest that 
the orbital period of the third component is longer than 1200- 
1500 d. Regular future eclipse timings will therefore be crucial 
in sorting out different mechanisms of period change and, ulti- 
mately, improving our understanding of the system. 
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Fig. 8. Radial velocity curve of DX Tuc. 



3.3. DX Tucanae 

The variability of DX Tuc was discovered by the Hipparcos 
satellite and it was classified as an F7-type contact binary (ESA 
1997, Kazarovets et al. 1999). Pribulla, Rreiner & Tremko 
(2003) listed the star in their catalog of field contact binary 
stars. Selam (2004), based on the Fourier-decomposition of the 
Hipparcos light curve, put it among the 64 genuine W UMa-type 
variables, while Pribulla & Rucinski (2006) did not find any in- 
dication for multiplicity (caused mainly by the lack of available 
data in the literature). 

The radial velocity curve and standard photometry presented 
here (Fig. [8][9j are the first ones in the literature. Using the 
Hipparcos ephemeris (HJD m ; n = 2448500.2540) and the new 
times of minimum light we calculated an improved orbital pe- 
riod P orb = 0.3771 1010(2) d. 

The full set of physical parameters is given in Table [5] The 
system is a typical contact binary, which has two components 
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QYHya 



Fig. 9. Observed and fitted light curves of DX Tuc. 



of very different masses but similar temperatures. Nevertheless, 
the larger star is slightly hotter, which puts DX Tuc among the 
A-type W UMa systems. The light curve shows a small but de- 
tectable O'Connell-efFect (AV = 0.01 mag), so we included a 
spot in the light curve model. With no indication of a third com- 
ponent in the cross-correlation profile, we fixed Z3 = 0. For cal- 
culating the absolute parameters, we assumed zero interstellar 
reddening, which is supported by the reddening map of Schlegel, 
Finkbeiner & Davis (1998) that implies E(B - V) < 0.017 mag 
in this direction. 
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3.4. QYHydrae 

The variability of QY Hydrae was discovered by the Hipparcos 
satellite (ESA 1997) and the variable star designation was given 
by Kazarovets et al. (1999). Although it belongs to the 100 
brightest X-ray stars within 50 parsecs of the Sun (Makarov 
2003), no detailed study has been done so far. (We note that 
Makarov (2003) listed the star erroneously as QY Lyr.) The late 
spectral type (K2V) and the short orbital period put QY Hya in 
the XO class of X-ray active stars by Makarov (2003), whose 
group includes binary stars of BY Dra-type, detached binaries 
(Algols) and eclipsing binaries of B Lyr type. Selam (2004) 
clearly separated QY Hya from the W UMa stars using the 
Fourier-description of the Hipparcos light curve, confirming the 
B Lyr class. There is one measurement of the radial velocity, 
+25.4 ± 0.6 km s~\ published by Nordstrom et al. (2004). They 
observed the star among ~ 14,000 F and G dwarfs, thus avoiding 
the detection of possible time variability of the radial velocity. 

We obtained standard BV light curves with full phase cover- 
age and an excellent coverage of the radial velocity curve (Figs. 
fTOlfTll i. The cross-correlation profile shows the presence of a 
third light at V3 ~ +50 km s -1 (top panel in Fig.fTOl). Using the 
Hipparcos ephemeris (HJD m j n = 2448500.2490) and our new 
epochs of minimum light, we determined an updated mean pe- 
riod of P orb = 0.29234050(8) d. 

The light curve fit indicates a semi-detached binary of simi- 
lar components (Tabled There is a slight but significant bright- 
ness difference of AV = 0.035 mag between the two maxima, 
which can be explained by a dark spot on the primary compo- 
nent. The calculated absolute parameters are based on the as- 
sumption of negligible interstellar reddening (the reddening map 
of Schlegel, Finkbeiner & Davis (1998) gives an upper limit of 
E(B - V) < 0.07 mag towards QY Hya, so the ~50 pc distance 
to the star indeed implies a small colour excess). The third light 



Fig. 10. Top panel: the CCF-profile of QY Hya in phase 0.75 
with the three Gaussian components. Bottom panel: radial ve- 
locities of the eclipsing components of QY Hya. 
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Fig. 11. Observed and fitted light curves of QY Hya. 



seems to have similar contributions in B and V, like the third 
star in VZ Lib, which suggests a hierarchic triple system of three 
similar components. 

3.5. V870Arae 

V870 Ara is another Hipparcos discovery, classified as an F8- 
type contact binary (ESA 1997, Kazarovets et al. 1999). Pribulla, 
Kreiner & Tremko (2003) listed the star in their catalog of 
field contact binary stars. Selam (2004), based on the Fourier- 
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Table 6. Physical parameters of QY Hya. Underlined values 
were calculated and fixed by the WD-code. HJD3 and V3 refer to 
mean epoch and radial velocity of the third component. 
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Parameter 


Value 


Spectroscopy 
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parameters 






Co-lat. [deg] 


33 


Rad. [deg] 


46 


Long, [deg] 


350 


Tfact 


0.84+0.01 
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U [L Q ] 
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Sp. type (3) 
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decomposition of the Hipparcos light curve, put it among the 
64 genuine W UMa-type variables, while Pribulla & Rucinski 
(2006) did not find any indication for multiplicity (caused largely 
by the lack of data in the literature). 

The radial velocity curve and standard photometry presented 
here (Fig. [T2"1[T3| > are the first ones in the literature. Using the 
Hipparcos ephemeris (HJD m ; n = 2448500.1840) and the new 
times of minimum light we calculated an improved orbital period 
p 0[h =0.39972200(2) d. 

The most interesting feature about V870 Ara is the very 
small mass-ratio q = 0.082 + 0.030. There are only two contact 
binaries with spectroscopically measured mass-ratios around or 
below 0.08: AW UMa (q « 0.075, Rucinski 1992), SX Crv 
(q ~ 0.066, Rucinski et al. 2001); and one star, V857 Her, for 
which the best-fit light curve solutions strongly suggest a mass- 
ratio somewhat less than 0.07 but lacking spectroscopic confir- 
mation (Qian et al. 2005). The existence of these stars are impor- 
tant, because theory currently predicts a cutoff at about q = 0.09 
(Rasio 1995), which might be pushed a bit lower to q — 0.076 
(Li & Zhang 2006). Below that contact binaries are expected 
to merge a single fast-rotating star within 10 3 -10 4 yrs. This puts 
V870 Ara among the objects that have the potential of constrain- 
ing evolutionary scenarios of binary mergers. 

Because of the well expressed O'Connell-effect (AV = 0.032 
mag) and asymmetric distortions of the light curve we added 
two spots to the light curve solution. They rather represent the 
difficulties we met during the light curve modelling than two 
real compact features on the hot component. The finally adopted 
set of parameters (Table [7]) shows that despite its extreme mass- 
ratio, this is a typical W-type W UMa system. 



Table 7. Physical parameters of V870 Ara. 
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Fig. 12. Radial velocity curve of V870 Ara. 

4. Summary 

In this paper we presented new photometric and spectroscopic 
data and their basic analysis for five close eclipsing binary stars. 
The sample consisted of three southern and two equatorial vari- 
ables, of which the southern objects have never been observed 
and modelled since the discovery. The main results of this inves- 
tigation can be summarized as follows: 

- XY Leo is a hierarchic quadruple system with a W UMa-type 
contact binary and a BY Dra-type red dwarf binary. It is one 
of the best cases for the light-time effect in a periodic vari- 
able star. Besides determining new spectroscopic elements 
and a light curve solution, we also found weak evidence for 
short-period magnetic cycles. 

- VZ Lib is another multiple system in which we detect the 
third component both spectroscopically and photometrically. 



10 



T. Szalai et al.: Physical parameters and multiplicity of five southern close eclipsing binaries 



9.2 - 



9.4 - 



9.8 





B 




V 











0.2 



0.4 



0.6 O.f 
Phase 



1.2 



Fig. 13. Observed and fitted light curves of V870 Ara. 
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Fig. 14. A comparison of the calculated distances with the 
Hipparcos measurements. 



Recent data indicate the possibility of detectable light-time 
effect, thus further eclipse timings are needed to measure the 
orbital period of the tertiary companion. 

- DX Tuc is a typical A-type W UMa star (i.e. the larger com- 
ponent is hotter). 

- QY Hya, being one of the 100 brightest X-ray active stars 
within 50 parsecs of the Sun, is a late-type triple system with 
a semi-detached eclipsing pair. 

- Finally, V870 Ara is a contact binary with the third smallest 
spectroscopic mass-ratio in all W UMa stars to date. 

The consistency of the presented results can be tested with 
the Hipparcos distances. A comparison of the calculated dis- 
tances that are based on the light curve models and the parallax- 
based Hipparcos values show a good agreement for all the pro- 
gramme stars (Fig.[T4l>. In case of DX Tuc and VZ Lib, the dis- 
tances from light curve models are likely to be improvements 
over the Hipparcos values. We note that both the two most 
deviant stars, XY Leo and VZ Lib, have bright tertiary com- 
ponents, which may have introduced a systematic error in the 
Hipparcos astrometry that can explain the larger disagreement. 
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